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ABSTRACT 

Elastohydrodynamic oil film measurements for rolling point con- 
tact under starvation conditions are obtained using optical interfer- 
ometry. The experimental measurements present a reasonably clear 
picture of the starvation phenomenon and. are shown to agree with 
theoretical predictions. Starvation inhibits the generation of 
pressure and, therefore, reduces film thickness. It also causes the 
overall pressure, stress and elastic deformation to become more 
Hertzian. 

Additional experiments using interferometry illustrate: , The 
.cavitation pattern, lubricant entrapment, grease lubrication, ball 
spin, and edge effects in line contact. 


INTRODUCTION 

The thickness of the hydrodynamic film between bearing surfaces 
is a very important consideration, particularly for elastohydrodynamic 
(EHD) conditions where film thickness and surface irregularities are 
often of the same order of magnitude. The present state of EHD shows 
that the development of an EHD film is now quite well understood, and 
that the thickness of such films, can be predicted theoretically with 
commendable accuracy. Presently, some investigations are concerned 
with isolated aspects of EHD and the explanations of discrepancies be- 
tween theory and experiment. Considerable emphasis is being placed on 
the application of EHD principles to establish the most favorable op- 
erating conditions in machine elements. 

Many applications, however, do not operate under the ideal con- 
ditions which are normally assumed or provided in most theoretical and 
experimental investigations. An important practical consideration 
which is often omitted is the supply and distribution of lubricant in 
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the vicinity of the contact region. Lubricant starvation, which may 
often exist but not always be recognized, can have an overriding 
effect on film thickness and other EHD aspects. 

An important aspect inherent in the EHD problem which was brought 
out in the Grubin theory [l] is that film thickness is primarily 
developed in the inlet region, just in advance of the flat Hertzian 
region. The experimental work of Crook [2] and the thermal solutions 
of Cheng and Sternlicht [3] and Dowson and Whitaker [4] provide fur- 
ther evidence of the fact that film thickness is governed by the con- 
ditions in the inlet region. Therefore, the conditions of speed, 
viscosity, and geometry, as far as film thickness is concerned, must 
b.e those relevant to the inlet region. These are the controlling 
variables in the film thickness equations. 

A condition of the inlet region which is not represented in the 
film thickness equations is that of the supply of lubricant. Hydro- 
dynamic considerations require the gap between the surfaces to b.e 
adequately filled to facilitate the important pressure buildup. If 
the inlet region is adequately filled with lubricant, the inlet 
pressure buildup will develop a film thickness which is insensitive 
to the further increases in lubricant supply. This condition of 
lubrication is defined as "flooded". If the inlet region is inade- 
quately filled, the inlet pressure buildup will be delayed, and a 
film thkckness which is dependent upon the available lubricant supply 
will result. The condition of lubrication is then defined as "starved". 
This starved condition of lubrication, which has received very little 
theoretical and experimental consideration, forms most of the work 
presented here. The experimental and theoretical work is limited to 
EHD point contact conditions; however, the fundamental features have 
more general application. 
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film thickness at. lubricant boundary 
Hertzian deformation for dry contact 
minimum film thickness 
central film thkckness 

central film thickness under flooded conditions 


dimensionless film thickness = h/h 

* ° 

transformed distance, m = log^ r 
pressure 

maximum Hertzian pressure 


reduced pressure =“(1-6°^) ? 

a h q 

o 

dimensionless reduced pressure = ^ 


maximum dimensionless reduced pressure = 
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radius from center of Hertzian region 
dimensionless radius = r/a 


111 

reduced radius of contact, — = + r~ 

Ja ^2 

radius of contacting bodies 

inlet distance, S = (x^ - a) 

inlet distance to obtain a flooded condition 

combined surface velocity = ^( u i + u p) 

velocity of bearing surfaces 

load on bearing contact 

coordinates 

distance from the inlet lubricant boundary to the center 
of the Hertzian region 

pressure viscosity coefficient 
angular dimension from axis of symmetry 
local viscosity of lubricant 
ambient viscosity of lubricant 
Poisson's ratio of bearing materials 
transformed angular dimension, = 20 
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Reynolds Equation 
Reynolds equation is 


THEORY 



Viscosity is assumed to vary with pressure according to 


A reduced pressure, q, is defined as 


(1) 


( 2 ) 


q = 


i 

a 


(1 - 


(3) 


Reynolds equation then becomes 


!r( h3 £) + M h3 iO ' 12u| vt (4) 

Grubin [l] applied this to line contact (where S'q/dy = o) by 
assuming (i) the shape outside the contact was Hertzian and (ii) 
that q « l/a everywhere inside the contact. The validity of this 
has previously been discussed (e.g., ref. 5)- The Grubin analysis has 
been extended to point contacts by Cameran and Gohar [6], Archard and 
Cowking [7], and more recently by Cheng [8]. 

The contact geometry is shown in Fig. 1. The film thickness h 
at any point is the sum of the central film thickness h and the 
Hertzian deformation for dry contact [ 5 ]. Thus 



(5) 
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Define 



To reduce the parameters of the problem the following terms axe defined: 
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Also, as the pressure falls rapidly with distance a new distance m 
is chosen m = log r*. Also, an angle 0 , chosen so it ranges from 
+1 to-1 is defined e as tt^ = 29 where 9 is the angular displacement 
from the line of fluid motion on the axis of symmetry. 


Applying these transforms Reynolds equation . takes the form 


cl 2 q* , 4 ^q 2 . 3 &h $q* 

5 m 2 * * * / ^ 2 h ^ »r 



1 &h 
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(8) 


Boundary Conditions 

1. Along AB (Fig. l) 

Due to symmetry along AB, = 0 


2 . Along BC 

q* = 0 along BC at some value of r*. This is the inlet boundary 
condition which is located at decreasing values of r* to obtain 
starvation conditions. 

3. Along CD (side leakage) 

q* = q* /h 2 along CD. The power of l/h was initially set to 

2 as in max [6]j but, for low values of B in the neighborhood of 

the Hertz radius minima would occur inside the boundary 0=1 in- 

stead of on the boundary. A more compatible boundary condition could 
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be obtained by increasing the power of l/h. However, the nature of 
the side leakage had no effect on the values along the axis of symmetry. 


4. Along DA (Hertz radius) ^ 

The Grubin solution shows q -4 - and d.q/dm 0 at the Hertzian 
radius. Therefore, q* = al$ng DA; however c4q*/^ m = 0 was 

only held at point A. The 19 ^* values along the axis of symmetry 
were used to estimate the derivative of q* at the Hertz radius 
(m = 0) . Sufficiently good starting values for q# at m = 0 re- 
sulted in derivatives which were linearly dependent on the q*^ 
value and thus linear interpolation of a few values enabled 
q*^ with dq*/8m = 0 to be attained. 

Finite - Difference Approximation 


The (m 5 <t>) plane is covered by a moderately fine mesh of up to 
50 equidistant steps in the m direction and up to 36 .in the <t> 
direction totaling 1800 points. 2 and 3 point finite differences were 
used as approximations to the partial derivatives yielding a 5-point 
difference equation for each mesh point. Application of such a 
formula to all mesh points results in a matrix equation Ax * b where 
x is the column vector of values at each mesh point, b is a column 
vector of constants and A is a (1800 x 1800) matrix. 


Inversion of a Large Matrix 

Such a matrix would be impossible to invert (or even store) on 
present-day computers. The nature of A is such that it may be 
partitioned into submatrices- w ith each row of submatrices consisting 
of three matrices (maximum 36 x 36 ), two of which are diagonal matrices. 

Inversion of the original matrix is achieved by Gaussian Elimina- 
tion Process with the submatrices and their inverses replacing the con- 
ventional elements and reciprocals. 

Starting Values 

A "reasonable" value for q^^ at the Hertzian radius (m = 0) 
for a given value of B is useatnd the m-interval (initially set so 
that the boundary is well outside the pressure field) is reduced thus 
bringing in the boundary until the vector of q- values beings to change. 
This then establishes the boundary for the particular flooded case 
and the minimum interval is achieved. Next the number of intervals, 
which is initially 36 , is reduced until the "picture" changes. The 
speed of the program is mainly determined by the size of this number. 
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Results for Flooded Condition 

The values of q* were determined for a fairly extensive 
range of B values wTtn log 1Q B ranging from -3 to +3. A plot of 
log, 0 <3* against log, „ B gave a good straight line in the range 
log 10 '°' 2 to 1 ”° Siting 


Snax " °*° 93 B ' 


• 0.6 


( 9 ) 


Remembering that at the Hertzian radius 



and 




12p o uaa 


( 10 ) 


the following film thickness formula for a flooded condition is 
obtained: 



Results for Starved Condition 

The results for a starved condition are represented in Pig. (2). 
Here q* I s plotted against n (where r* = 1 + 2 n "' ) ) for con- 
stant values of q*B . q*B can be determined from the operating condi- 
tions as shown in [ 9 ]. - 



For a given inlet condition r*, q* can be determined from Fig. (2) 
and h Q can be calculated from equation (10). 
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EXPERIMENTAL PRELIMINARIES 


Apparatus 

A schematic diagram of the apparatus is shown in Fig. 3* The 
apparatus provides a system whereby the EHD film thickness in a ball 
thrust bearing can be measured using optical interferometry. The 
experimental bearing is shown in Fig. 4. It consists of two flat 
races separated by three 1-inch diameter steel balls. The lower race 
is made of steel and is supported by an air thrust bearing which is 
driven by an electric motor. The upper race is an optical crown 
glass disk 6 inches in diameter and l/2-inch thick. Its lower surface 
is coated with a thin (~170 X) chromium layer of 15 percent reflectivity, 
and its upper surface has an antireflection coating. Load is applied 
to the upper race with a lever arm, and the interference fringes are 
observed with a microscope having a magnification of X 50 . 

The bearing was operated in a counterrotation mode for most of 
the experiments by holding the cage stationary (Fig. 4). This causes 
the test ball, which can be adjusted radially in a laminated phenolic 
cage, to rotate about a single axis fixed in space. Both bottom and 
top race speeds were measured with a magnetic transducer. Slip 
measurements between races were always less than 2 percent for all ex- 
perimental conditions. 

The bearing was also operated by holding the top race stationary. 
The balls then have curvilinear as well as rotational motion. The 
interference fringes of the test ball are then observed stroboscopically 
by synchronizing the light source with the cage by means of a lamp and 
phototransistor as shown in Fig. 3- 

The microscope is illuminated with a xenon discharge tube which 
emits white light at a color temperature of 6000 K. The tube can be 
discharged at a high repetition rate for continuous observations, and 
it can be discharged singly for high-speed photography. 

The requirements for good fringe "visibility" and how they can be 
used to design an interferometric system which is compatible with the 
experimental measurements is discussed in [ 9 ]. A duochromatic system 
which is composed of green and red wavebands, obtained with a Wratten 
77A filter, provides fringes of good "visibility" over a satisfactory 
film thickness range of 80 microinches. The fringes are closely spaced 
and form a definite color sequence which can be used to determine 
fringe order. The calibration of this system is given in [ 9 ] and [10]. 

To obtain the absolute thickness corresponding to each fringe the 
calibrated fringes (obtained in air) must be divided by the refractive 
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index of the medium. A maximum experimental pressure of 87 000 psi 
increases the refractive index of a typical oil by about 8 percent 
from its atmospheric value. The Lorenz-Lorentz relation [11] was 
used to determine the change of refractive index with density, and 
the Hartung empirical formula for hydrocarbons [12] was used to deter- 
mine the change in density with pressure. 

Film Thickness Measurements (Flooded) 

Experimental film thickness measurements using a lubricant (LUB D) 
that is well documented [13 ] are plotted in Fig. 5 a function of 
four theoretical film thickness formulas. There is a very good agree- 
ment between theory and experiment for the present formula (Eq. 11) 
and that of Cheng [8]. These solutions were based on Grubin assump- 
tions have enabled a reasonably accurate film thickness formula for 
point contacts to be developed in the same manner as they did for 
line contact in the 1950 f s when a full solution had not yet been 
achieved. 

Figure 6 is a plot of the measured film shape obtained from a 

high-speed photomicrograph [9]. The measurements were corrected for 

refractive index change by using an estimated pressure profile. The 

film shape in the transverse direction shows that the minimum film 

thickness h at the lateral constrictions (h =8.5 microinch) is 

much less th§n the central film thickness (h = 18.3 microinch). The 

ability to predict the minimum film thickness is of practical importance. 

The central film thickness (h ) and the minimum film thickness (h ) 

' o HI 

as measured from random photomicrographs taken over a range of loads 

between I.92 lbf and 12 lbf are plotted in Fig. 7- The data for a 

materials combination of sapphire/tungsten carbide given by Gohar 

[14] are also included. The plot shows a certain degree of scatter 

because h is more sensitive to load than h . However , to a first 

approximation 


h 


m 


h -4.5 microinches 
o 


'(13) 


An interesting feature of the results is that the load parameter 
for the sapphire/tungsten carbide is nearly the same as the 
lowest load parameter for the glass/steel data. These two cases show 
very good correlation even though the reduced elastic modulus of the 
former (66.6x10° psi) is four times the latter, and the maximum Hertzian 
pressure of the former (405 000 psi) is over eight times the latter. 
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STARVATION EXPERIMENTS 
Inlet Lubricant Boundary and Film Shape 

The starvation experiments were conducted with the bearing oper- 
ating in the counterrotation mode. Before each test the bearing com- 
ponents were cleaned with a solvent (either benzene or acetone). 
Lubricant was applied to the ball and raceway surfaces during the re- 
assembly of the bearing. Under static conditions surface tension 
causes the lubricant to form a circular boundary or meniscus around 
the contact region. A side view of the meniscus can be seen in Pig. 4. 

Under dynamic conditions the quantity and distribution of lub- 
ricant in the inlet region is reflected in the share and location of 
the lubricant boundary in the inlet region. Figs. 8(a) and (b) show 
clearly defined areas where the gap between the bearing surfaces is 
completely filled with lubricant. Since the medium is continuous in 
these areas they form very uniform interference fringes. Outside 
these areas the medium is a composite of air and oil, and the fringes 
which form between the ball and glass race are shifted to locations 
of greater thickness because the refractive index of air is less than 
oil. This is clearly seen in the inlet region of Fig. 8(b). The 
degree of fringe shift depends on the relative amounts of air and oil 
at each location. In this way the rippled fringes in the inlet 
region show quantitatively the distribution of lubricant on the bear- 
ing surfaces. A fringe which is shifted closer to the Hertzian region 
indicates a greater percent of oil. 

The rippled fringes in the inlet region also show that the lub- 
ricant is distributed over the bearing surfaces in the form of small 
ribs which are generally oriented in the direction of motion. These 
ribs are developed in the outlet region where the lubricant film 
ruptures into striations of air and oil. These ribs, though somewhat 
flattened out, can also be seen in the inlet region. In this way 
the amount and distribution of lubricant in the outlet region influences 
the inlet region. 

Excess lubricant is frequently found on the sides of the track 
as shown in Fig. 8(b). These areas act as reservoirs from which lub- 
ricant can flow into the inlet region by the action of surface tension. 
The rate of flow is probably a function of surface tension, viscosity, 
gap thickness, and the shape of the lubricant boundary; and the amount 
of lubricant which is recaptured in the inlet region by surface ten- 
sion will certainly depend on the time available, i.e., the speed of 
the bearing surfaces. The amount and distribution of lubricant which 
fills the gap at any given time is clearly determined by the inlet 
lubricant boundary. 
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The EHD films shown in Fig. 8 are clearly affected by a starved 
inlet region^ and the film shape, which can be obtained by reading 
the fringes as lines of constant thickness, is strongly dependent 
upon the shape of the inlet lubricant boundary. The EHD film of 
Fig, 8(c) is very severely starved and consists of a series of 
channels. These channels are caused by the ribs of lubricant which 
can be seen in the inlet region. 

These results not only show that the film thickness is primarily 
determined in the inlet region, but also that the local film thick- 
ness at each point within the Hertzian region is primarily a function 
of the micro inlet conditions immediately upstream of that point. 

Test Conditions and Procedure 

The purpose of the starvation experiments was to determine the 
effect of the inlet boundary on the central film thickness for a 
range of operating conditions. For constant speed, load, and lubricant 
properties the required information is the central film thickness h 
and the boundary distance x^ measured from the center of the ° 

Hertzian region to the inlet boundary immediately upstream (see Fig. 9)* 
These measurements were taken from high-speed photomicrographs. 

Starvation' was achieved by gradually reducing the quantity of 
lubricant within the bearing system. The chosen operating conditions 
consisted of three loads on the test ball (1.92, 6, and 12 lbf) and 
four initial (i.e., flooded) film thicknesses which correspond to the 
first four red fringes of the duochromatic system. All combinations 
of the chosen operating conditions were tested, thus giving twelve in- 
dividual test runs. The three test loads give a maximum Hertzian 
pressure ranging from 47 400 to 87 000 psi, and the four initial film 
thicknesses range from 8 to 30 microinches. To maintain a convenient 
operating speed, on the order of 10 in. /sec, these film thicknesses 
were obtained by varying the test lubricant as well as the bearing 
speed. Five test lubricants were used. Their properties are shown 
in Table 1. The lubricants designated SAE 10, SAE 30, and SAE 40 
are paraffinic mineral oils, while TN 631 is a medium viscosity index 
naphthenic mineral oil. LUB D is a high viscosity, low viscosity 
index mineral oil which was previously used by Dyson et al. [15] and 
Foord et al. [16]. It was used here as a reference fluid from which 
the pressure-viscosity coefficients of the other test lubricants were 
obtained. The apparatus itself was used to obtain the pressure- 
viscosity coefficients in the manner described by Foord et al, [16]. 
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Starvation Results 


Figure 10 shows a sample of photomicrographs taken from a test 
run. The test load is 1.92 lbf, and the initial film thickness corres- 
ponds to the fourth red fringe (h =30 microinches). In Fig. 10(a) 
the boundary distance is sufficiently removed from the center of the 
Hertzian region to obtain a flooded condition. Fig. lO(b-f) show a 
diminishing film thickness as the inlet boundary progressively 
approaches the Hertzian region. The diminishing film thickness is 
apparent if one observes that the fringe in the center of the Hertzian 
region for a given figure is always in the constriction regions in 
the previous figure. It should also be noted that as starvation pro- 
gresses , the EHD shape approaches the Hertzian shape. Therefore, the 
EHD pressure must also approach the Hertzian pressure. In the limit, 
i.e., when h equals zero, the Hertzian shape and pressure are 
attained. 

Some of the experimental and theoretical central film thicknesses 
are plotted in Figs. 11 and 12 as a function of x,/a where a is 
the calculated Hertzian radius. Figure 11 shows the effect of various 
initial film thicknesses for a given load, and Fig. 12 shows the 
effect of load for a given initial film thickness. The results in- 
dicate that film thickness becomes increasingly more sensitive to the 
inlet boundary as the boundary approaches the Hertzian radius. Indeed, 
the film thickness approaches zero as the inlet boundary approaches 
the Hertzian radius. This emphasizes the importance of the inlet 
region in regards to film thickness. To obtain a flooded condition 
Fig. 11 shows that the required inlet boundary, expressed as the ratio 
x^/a, increases with film thickness. This reflects the expected 
inlet pressure buildup which commences further upstream as the film 
thickness (and hence speed and/or viscosity) increases. Figure 12 
shows the required inlet boundary x,/a to decrease with load. This 
reflects the overall pressure distribution which must become more 
Hertzian as the load increases. 

The agreement between theory and experiment is quite good when 
one considers the scatter in the experimental data together with the 
Grubin assumptions inherent in the theoretical solution, It must 
also be remembered that the film thickness predicted by theory is sub- 
jected to the accuracy of the measured experimental operating con- 
ditions, particularly the lubricant properties of a and p . It is 
of interest that the agreement between theory and experiment is no 
worse for starvation conditions than it is for flooded conditions. 

In fact, it is likely to be better under starvation conditions since 
the film shape becomes more like the assumed Hertzian shape when the 
inlet region is starved. 
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Prediction of the Onset of Starvation 

The above results imply that film thickness is determined by 
the EHD pressure generated in the inlet region and that starvation 
will initiate when the inlet lubricant boundary approaches the ex- 
pected location of pressure commencement. The approximate location 
of pressure commencement can be obtained from the integrated form 
of Reynolds equation. Neglecting the side~leak.age term equation 
(4) can be integrated to give the reduced pressure gradient 



For constant speed and viscosity the reduced pressure gradient with- 
in the inlet region is determined by the geometry term (h - h )/h^. 
This term, and therefore the reduced pressure gradient, reaches a 
maximum when h/h =1.5 and diminishes approaching zero in the up- 
stream direction as h/h increases. If h^ is the thickness of 
the lubricant film at the inlet boundary (see Fig. 9)> then it is 
expected that large values of the inlet geometry ratio li/h will 
provide a flooded condition. This is shown in Fig. 13 where all the 
experimental data are plotted as h /( h )„ against h^/h where 
(h ) is the film thickness under flooded conditions. The solid 
curve represents the theoretical results over the experimental range 
of operating conditions and shows very good agreement with experiment 
except near the flooded condition. The interesting feature of the 
results is that the degree of starvation, over a range of loads and 
initial film thicknesses, can be described by a single curve. The 
results also show that the flooded condition is approached asymptot- 
ically as increases. Therefore, the onset of starvation 

cannot be very well defined. However, most of the experimental re- 
sults indicate a flooded condition at h^/h = 9> and the theoretical 
curve predicts a film thickness which is 95 percent flooded. 

For practical applications a more appropriate measure of the 
degree of starvation is the inlet distance S defined as 


S = (Xh - a) 


(15) 


and shown diagramatically in Fig. 9. S can be written in terms of 

h../h as follows: By definition 

b o 
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\ h o * ^ 
h h 


(16) 


where Ir. is the Hertzian deformation for dry contact (Fig. l) . 
For poirS contacts can be approximated by (see ref. 5) 



This relation is accurate only for small ratios of x^/a. Figure 11 
shows that over the experimental range starvation occurs within 

x,/a =4.5. The calculation of from equation (17) for 
x!Va =4.5 is about 14 percent less than that calculated from the 
Hertzian relation (Eq. 5). For most EHD conditions starvation occurs 
within an x_/a ratio of 3 or 2 where the error is 7.6 percent and 
2 percent, respectively. Now we can write 


la / a 

and 


p 

max _ _a 
E' ~ rR 

Substituting equations (17), (18) and (19) in equation (16) and solv- 
ing for S gives 


(18) 


(19) 



(20) 
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A similar relation can be derived for EHD line contacts [9]. If 
the onset of starvation is defined by h^/h = 9 s then the minimum 
inlet distance to obtain a flooded conSition can be approxi- 

mated by 


3„52(R h ) 2//3 

Q - ° 

S f a l/3 


( 21 ) 


This is very similar to the statement of Fein and Kreutz [17 ] that 
the vast majority of pressure buildup for an EHD point contact occurs 
within about 


3 (R h o ) 

~“T7T 


2/3 


ahead of the Hertzian flat. 

Figure l4 is a plot of all the starvation results using the in- 
let distance S as a measure of the degree of starvation. The ex- 
perimental data is represented by four curves each showing the film 
thickness approaching zero as the inlet distance approaches zero. 

Each curve contains the data for all three test loads thus indicating 
that starvation is not very sensitive to load. This is consistent 
with the view that film thickness is primarily determined by the EHD 
pressure generated in the inlet region. The primary effect of load 
is merely to alter the size of the Hertzian region and the magnitude 
of the pressure therein. Load has only a secondary influence on the 
EHD pressure generated in the inlet region. Therefore;, this pressure, 
and hence the film thickness determined by this pressure, will re- 
main essentially the same regardless of load. 

The onset of starvation given by equation (21) is plotted in 
Fig, l4 for each test load. It compares very well with the experi- 
mental results. An interesting feature is that for a given central 
film thickness h , equation (21) predicts that starvation will com- 
mence at slightly°smaller inlet distances when the load increases. 

The experimental data also reveals this tendency. This variation with 
load occurs because the gap thickness increases with load for a con- 
stant S and h , thus making the ratio h^/h^ larger. 

Also plotted in Fig. l4 is the line defining the equation 


S = 500 h Q 


( 22 ) 
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Over the experimental operating conditions this equation approximates 
the onset of starvation and indicates the relative magnitude of the 
dimensions involved. 

Semi -Empirical Formula for Starvation 

All the starvation results are replotted in Fig. 15 where s/s f , 
as a dimensionless starvation parameter, now represents the degree 
of starvation. The experimental data can be approximated by a quadrant of 
a circle, i.e.. 



This can be employed in equation (ll) to provide the following semi- 
empirical film thickness formula which includes the effects of starvation 



The greatest practical problem in applying equation (24) is the 
determination of the inlet distance for the particular application. 

The inlet distance is difficult to measure directly and difficult to 
predict from external conditions. Besides the method and amount of 
lubricant supplied, these external conditions must involve those aspects 
which deal with the transport and distribution of lubricant within the 
bearing, i.e. viscosity, surface tension, bearing dynamics, and geometry. 
There is certainly a great need for practical research along these lines. 


EFFECT OF STARVATION ON OTHER ASPECTS OF EHD 

Friction 

Friction arises from the viscous shear forces acting on the bear- 
ing surfaces. For EHD contacts rolling friction is associated with the 
pressure generation in the inlet region, and sliding friction, or traction, 
is associated with the flat Hertzian region and is inversely proportional 
to the central film thickness h . A starved inlet region inhibits the inlet 
pressure buildup and reduces the°centra! film thickness. Therefore, a 
reduction in rolling fhiction and an increase in sliding friction, or 
traction, is expected to accompany starvation. Both these effects have 
practical implications. In rolling element bearings, for example, it is 
desirable to have low rolling friction and high traction. The latter is 
important for preventing slip, and hence bearing failure, between the 
rolling elements and raceways. 
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Pressure and Stress Distribution 

It was shown in Fig. 10 that as starvation progresses the EHD 
film shape approaches the Hertzian shape. Thus, it should follow that 
the pressure distribution and therefore the stress distribution will 
also approach the Hertzian condition. 


ADDITIONAL EXPERIMENTS 


Cavitation 

Ordinary lubricants contain about 10 percent (by volume) of dis- 
solved gas or air at atmospheric pressure. When the lubricant pres- 
sure is reduced below the saturation pressure in the exit region the 
dissolved gas diffuses into small gas bubbles within the lubricant 
causing the bubbles to grow in size. This form of bubble growth, 
described by Schweitzer and Szebehely [18] and Floberg [19J } is the 
usual form of cavitation which occurs in lubricating films. 

The cavitated region, under slow speeds, may often be contained 
within the lubricant meniscus surrounding the contact region. At 
very slow speeds it may consist of a single gas bubble like that ob- 
served by Ax chard and Kirk [20]; or, at higher speeds it may develop 
into a large number of gas bubbles. If the bearing is stopped, the 
bubbles often, coalesce into a single bubble which detaches itself 
from the contact region and slowly goes back into solution as it 
moves downstream. This is shown in Fig. 16. Figure l6( a), which 
was taken about one second after the bearing had stopped, shows a 
bubble just prior to its separation from the contact region. Fig- 
ure 16 (b) shows a bubble approximately two seconds after the bearing 
had stopped. An interesting feature is the interference pattern 
which has developed, apparently between the glass race and the sur- 
face of the gas bubble. It indicates an "EHD” film which has the 
features of the classical EHD point contact shape. Both the lateral 
and rear constrictions are visible as well as the slight increase in 
film thickness must prior to the rear constriction. 

Under normal bearing speeds the cavitated region breaks through 
the lubricant meniscus on the downstream side, thus exposing itself 
to the atmosphere. This has little effect on the general appearance 
of the cavitated pattern. The lubricant which passes through the 
Hertzian region generally ruptures at a film thickness which is a- 
bout 1.5 times the central film thickness. The cavitation pattern 
consists of thin ribs of lubricant separated by air spaces. This is 
most clearly seen for large central film thicknesses like that shown 
in Fig. 17(a) . 
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The lubricant which escapes around the Hertzian region ruptures 
into cavitation fingers which grow perpendicular to a line which 
makes an angle of 35 with the direction of motion. The shape of 
the cavitation fingers is partially determined by surface tension. 

A high surface tension fluid causes the cavitation fingers to become 
more rounded. This is shown in Fig. 17(b) where the lubricant is 
glycerol. It has a surface tension of 63 dynes/cm which is twice 
that of most lubricating oils. The small bubble to the right of the 
Hertzian region has cavitated prematurely. This may be due to a 
foreign particle within the fluid which has provided a preferential 
nucleating site for cavitation to take place. 

The lubricant in Fig. 17(c) is a concentrated dispersion of 
very fine colloidal graphite (Oildag) in SAE 30 oil. These graphite 
particles provi.de additional nuclearing sites thus causing the 
branches to form on the cavitation fingers. 

In most bearing applications a cavitated outlet region is 
immediately followed by an inlet region. The influence of the cavitated 
region on the lubricant supply condition in the inlet region is an 
important practical aspect which is often overlooked. It is of in- 
terest, therefore, to consider the distribution of lubricant within 
the cavitated region. 

The relative shift of the interference fringes due to the different 
refractive indices of air and oil indicate the distribution of lubri- 
cant within the cavitated region. Figure 18(a) shows schematically the 
distribution of lubricant within the cavitated region at a constant 
radius from the ball center (i.e., where the gap thickness is constant). 
The individual air spaces, as well as the ribs of lubricant, within the 
ball track essentially extend across the gap between the bearing sur- 
faces, After the ball has passed, the glass race is left with the 
situation shown schematically in Fig, 18 (b), Narrow ribs of lubricant 
remain on the track while thin waves of lubricant form adjacent to the 
track. The same general lubricant distribution, though somewhat 
flattened out, appears in the inlet region. Unless this track is filled 
in, it is this depleted supply of lubricant which must provide a 
sufficient quantity of lubricant in the inlet region to maintain a 
flooded condition. If the design of a bearing and the method of lub- 
ricant supply do not fulfill the local lubricant requirements in the 
inlet region, the bearing may be operating under starved conditions even 
though it is "well lubricated". 

Entrapment 

An oil entrapment between the ball and glass race is shown in Fig. 
19(a). The surfaces are in contact (or very nearly so) only along the 
edge of the Hertzian region. The fringes, which are of a monochromatic 
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green system, indicate a wedge-shaped film shown schematically in 
Fig. 19(b). 

The particular shape of the entrapment and its formation is some- 
what different from the normal approach entrapments described by other 
workers [16, 21-23]. The bearing was rotated slowly in a viscous lub- 
ricant (LUB D) and then stopped by applying a breaking torque. The 
wedge-shaped film is probably the result of a decreasing surface speed 
which generated a thinner film in the inlet region as the motion was 
stopped. The braking torque produced a slight reverse motion equal 
to about one-fourth the Hertzian diameter. This caused some of the 
fluid to leak out and produced the crescent-shaped ridge in the film 
thickness along the right-hand edge of the Hertzian region. 

The photomicrograph was taken approximately five minutes after 
the entrapment was formed. During that time some of the lubricant 
leaked out around the edge of the Hertzian region and along the 
scratches in the bearing surfaces. These results imply that in a 
bearing application where there is a velocity reversal it is possible 
to maintain a full EHD film even though the relative surface velocity 
is momentarily zero. 

Grease Lubrication 


Figure 20 is a plot of film thickness against speed for three 
lithium hydroxystearate soap greases and their base oils. The base 
oil viscosity and soap content of the greases are given in table 2. 
Each grease gives a higher film thickness than its base oil thus in- 
dicating that the apparent viscosity of the grease is greater than 
the base oil viscosity, even at very high rates of shear (10 sec -1 ). 

It is also noted that the grease G2 with 10.5 percent soap content 
gives a greater film thickness than the grease G3 with 6 percent soap 
content. This is in line with recent work by Dyson and Wilson [24] 
on the same greases which indicated that the soap breaks down into 
hard spheroidal particles. These particles have a thickening effect 
on the base oil which prevents the apparent discosity from reaching 
the base oil viscosity. 

A comparison of the grease and oil film shapes under moderate 
speeds showed no noticeable difference, provided they were compared 
at the same nominal film thickness. At extremely slow speeds, how- 
ever, relatively thick and irregular films were observed with grease 
particularly when operating in fresh grease. These films remained un- 
der static conditions as shown in Fig. 21. Figure 21(a) was obtained 
after the ball had been rolled very slowly in fresh grease (Gl). The 
fringes within the Hertzian region are very irregular in. shape, thus 
indicating a nonuniform film shape or possibly large changes in re- 
fractive index due to the heterogeneous structure of the grease. 
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Figure 21(b) was obtained after the ball had been rolled several times 
over the same track. The resulting film is thinner (note the position 
of the black fringe in each photomicrograph) and more uniform, thus 
suggesting a breakdown in soap structure. It is possible that the 
small circular spots in Fig. 21(b) are broken down soap particles. It 
was only possible to obtain a film of the thickness shown in Fig. 21(a) 
by resupplying the bearing with unworked grease. 

A test, run at constant speed, showed no change in film thickness 
with grease G1 over a period of 30 minutes. It is likely, therefore, 
that most of the soap structure had broken down soon after the bearing 
was started. The test was terminated after 45 minutes when the film 
thickness became clearly affected by starvation. This was always 
found difficult to avoid. 

Starvation tests performed with grease G2 and its base oil showed 
no apparent difference in the way the film thickness diminishes with 
the inlet distance from its initially flooded value. Important differ- 
ences, however, between grease^ and oil lubrication in connection with 
starvation are likely to be found in the mechanisms by which the 
lubricant is transported and redistributed in the vicinity of the in- 
let region. 

Ball Spin 


Figure 22 shows two photomicrographs taken with the top race of 
the thrust bearing held stationary. The quality of the image is lack- 
ing because at high object speeds the flash duration is not short 
enough to prevent image blur. In Fig. 22(a), where the cage speed is 
263 rpm, the film shape and wake pattern are similar to that found 
when the bearing is operating in a counterrotation manner. Figure 
22(b) was taken when the cage speed was increased to 300 rpm. Both 
the film shape and wake pattern are skewed on the order of 15° away 
from the orbital direction of the ball. This can be attributed to a 
gyroscopic moment which produces a sliding velocity component in the 
transverse direction. The gyroscopic moment in Fig. 22(a) is resisted 
by the frictional forces at the cage and bearing races. This re- 
quires an effective coefficient of friction of 0.02 for the particular 
operating conditions [ 9 ]* 

Edge Effects in Line Contact 

Preliminary work was carried out on a tapered roller thrust bear- 
ing where the top glass race was held stationary and the interference 
fringes were observed stroboscopically. In Fig. 23(a) the deformation 
of the glass beyond the roller end causes an additional compressive 
stress at the roller end which produces the characteristic "dog bone" 
shape. Also, a very predominant lateral constriction in film thickness 
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can be seen at the contact edge. This can be attributed to side- 
leakage effects, primarily in the inlet region. 

The edge effects can be reduced by blending the roller ends, as 
shown by Gohar and Cameron [ 25 ], or crowning the roller profile. A 
roller with a blended end is shown in Fig. 23(b). The central film 
thickness is the same as that in Fig. 23(a) but the lateral con- 
striction has almost been eliminated. 

Some interesting implications arise from the film shapes in 
Fig. 23. The film shape at the edge of the blended roller is very 
similar to that found in point contacts. Even the cavitation pattern 
along the edge is similar to the point contact case. In this respect, 
point contacts can be considered as " naturally" blended or crowned. 
While point contacts may not have the optimum crown to eliminate the 
lateral constrictions, it is possible that an elliptical contact found, 
for example, between a ball and grooved raceway may provide the re- 
quired geometry for its elimination. These results further imply 
that the lateral constriction for an unblended roller, as found in 
Fig. 23(a), may be far more severe than those associated with point 
contacts. 


CONCLUSIONS 

Starvation under EHD conditions is associated with the location 
of the lubricant boundary in the inlet region. Under starvation 
conditions the local film thickness in the Hertzian region is a func- 
tion of the lubricant ' boundary immediately upstream. The results 
show that the onset of starvation is reflected in the expected loca- 
tion of pressure commencement, and that film thickness diminishes to 
zero as the lubricant boundary approaches the Hertzian region. This 
implies that film thickness is established by the EHD pressure gen- 
erated in the inlet region. The primary effect of starvation is that 
it inhibits the generation of this pressure. The effect of this was 
measured in terms of film thickness and was found to be in harmony 
with EHD theory based on Grubin assumptions. 

Under constant operating conditions the hydrodynamic pressure 
generation is a function of the inlet gap geometry. The ratio h/h o 
describes the relative thickness of this gap and reflects the 
relative magnitude of the hydrodynamic pressure generation along the 
inlet region. The experimental results correlate with the inlet 
geometry ratio h^/h^ where h^ is the gap thickness at the lubricant 

boundary. The onset of starvation is not very well defined, but can 
be approximated by h b /h o = 9 . This can be used to derive the onset 
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of starvation in terms of the inlet distance S which is the dis- 
tance between the lubricant boundary and the edge of the Hertzian 
region. If the location of the inlet boundary is known, film thick- 
ness can be predicted theoretically or by a semi-empirical formula 
based on the experimental results. 

Film thickness is only one aspect of starvation. The suppression 
of hydrodynamic pressure by starvation causes the overall pressure 
and stress distribution to conform to the Hertzian condition. It 
also reduces rolling friction and increases traction. Because these 
aspects have practical implications, starvation should not be con- 
sidered as just a situation to avoid, but should be dealt with as 
an additional EHD parameter which can be usefully employed. 
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TABLE I. - PROPERTIES OF TEST LUBRICANTS 


Lubricant 

symbol 

Specific 

gravity 

Viscosity 

5 CS : 

Pressure-viscosity 

coefficient 

(psi" 1 ) 



210° F 

SAE 10 

0.870 

66.8 

33.9 

5.2 

1. 41x10" ^ 

SAE 30 

.883 

261.0 

113.0 

11.4 

1.60 

SAE 40 

.885 

401.0 

I 65 .O 

14.8 

1.70 

TN 631 

• 901 


138.0 

10.2 

2.14 

LUB D 

.942 

572.0 

182 

, 1, , ,1 

10.8 

2.50 


TABLE II. - PROPERTIES OF LITHIUM 
HYDROXYSTEARATE SOAP GREASES 

























Figure 1. - Contact geometry. 
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Figure 2. -Theoretical results for starvation. 







Figure 4. - Experimental ball thrust bearing with stationary cage. 
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SYMBOL SOURCE THEORETICAL CONSTANTS 
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° ARCHARD & COWKING 
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» CAMERON & GOHAR 
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Figure 5. - Comparison of point contact theories with 
experiment 
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Figure 6. - Measured film shape (a) along the centerline in the 
direction of motion, and (b) in the transverse direction at the 
location of minimum film thickness; u = 2. 74 in. /sec; 

U 0 = 80x10^ (lbf)(sec)/in. 2, a = 2. 5xl0' 4 psi'l, w ■ 12 Ibf, 
Pmax = 87 000 psi. 
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Figure 7. - Plot of the central film thickness versus 
minimum film thickness for a range of experimental 
conditions. 




Figure 9. - Contact geometry showing required meas 
urements for starvation; (a) plan view, (b) cross- 
sectional view along centerline. 


Figure 8. - Effect of inlet lubricant boundary on 
film shape. 
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Figure 10. - Photomicrographs showing progressive starvation 


Figure 12. - Effect of inlet boundary on film thickness for dif- 
ferent loads. 
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Figure 13. - Comparison of theory and experiment for all the 
starvation results using the inlet geometry ratio h b /h Q . 



Figure 14. - Experimental results plotted as film thickness 
against inlet distance. 
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Figure 15. - The degree of starvation as a 
function of the starvation parameter S/Sj. 
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Figure 16. - Coalescence and separation of the cavitation bubbles: 
(a) about one second after bearing stopped, (b) about two seconds 
after bearing stopped showing the "EHD lubrication" of an air 
bubble. 


Figure 17. - Cavitation pattern in rolling point 
contact: (a) for large film thickness (h Q ■ 67 


microinches) showing striations of air and oil, 
(b) for high surface tension fluid (glycerol), 
and (c) for a concentrated dispersion of 


colloidal graphite (Oildag). 




Figure 18. - Schematic representation of the cavitated region: (a) at 
a constant radius from the ball center, and (b) showing the distri- 
bution of lubricant after the ball has passed. 



Figure 19. - (a) Photomicrograph of an oil entrapment; (b) schematic 
representation of film thickness. 
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Figure 20. - Film thickness versus speed for three greases 
and their base oils. 


igu re 21 . - Grease films under static conditions: (a) ob- 
tained after rolling in fresh grease, lb) obtained after rol 
ing several times over the same track. 







Figure 22. - Observation of ball spin: (a) no 
spin at cage speed of 263 rpm, (b) ball spins 
at cage speed of 300 rpm. 



Figure 23. - Edge effects in line contact: (a) un- 
blended roller end, (b> blended roller end. 
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